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■ 

The  research  performed  under  this  contract  during  the  period  1  Oct.  1978 
through  31  March  1979  can  be  divided  into  three  main  topics,  upper  mantle 
studies,  propagational  distortions  caused  by  non-planar  structure  and  S  wave 
generation  by  block  movements  near  explosions. 

■2- 

In  Section  leT,  we  discuss  the  upper  mantle  structure  of  northwestern  Eurasia 
where  long  period  and  short  period  WWSS  seismograms  from  nuclear  explosions  in 
the  Soviet  Union  have  been  incorporated  with  published  apparent  velocity  data  — 

DD  i  j*"M73  1473  hmt.i  a'A'vtftfp 
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k-to  derive  a  comp r ess  Iona  1  wave  velocity  model  lor  the' upper  mantle  beneath  north¬ 
west  Kurasia.  The  waveforms  were  modeled  using  generalized  rav  theory  and  a 
starting  model  was  perturbed  by  trial  and  error  until  a  lit  was  achieved.  The 
seismograms  from  9-13°  exhibit  impulsive  first  arrivals  (Pn) ,  implying  a  smooth, 
positive  velocity  gradient  between  depths  of  60  and  150  km.  A  consistent  puls« 
arriving  about  2  seconds  after  P„  at  the  ranges  11  to  17“  is  modeled  by 
positioning  a  low  velocity  zone  between  150  and  200  km.  The  Pn  phase  becomes 
diffracted  beyond  15“  although  there  is  considerable  evidence  of  lateral 
variation  at  where  this  geometric  cut-off  occurs.  The  model  is  relatively 
smooth  from  a  depth  of  200  km  down  to  420  km  where  a  5%  jump  in  velocity 
produces  a  triplication  from  15  to  23°.  The  observations  from  21  to  26° 
clearly  show  another  discontinuity  at  a  depth  of  675  km  with  a  4%  change  in 
veloc l ty 

In  Section  3,  glorified  optics  has  been  used  to  compare  the  waveform 
complications  caused  by  two  and  three  dimensional  structures.  A  scheme  is 
developed  to  synthesize  seismic  images,  in  which  the  direct  rays  give  a 
preliminary  feature  of  a  structure,  while  the  multiple  boundes  provide 
further  information  for  Improving  and  refining  the  model.  A  method  is 
discussed  for  inverting  non-strati fied  structures  like  sedimentary  basins  and 
subduct  ion  plates. 

*  In  Section  4,  we  consider  the  problem  of  reflection  and  refraction  of 
purely  compress iona 1  waves  incident  on  an  interface  separating  identical  solid 
half-spaces  In  which  the  condition  of  continuity  of  shear  displacement  at  the 
boundary  is  generalized  to  one  that  allows  slippage.  The  problem  is  solved 
using  tile  Cagniard-de  Hoop  technique.  It  is  Found  /^hat  the  generation  of 
reflected  1’  and  S  waves,  as  well  as  transmitted  S  wnvVs,  is  most  effective 
In  the  case  of  a  perfectly  unbounded  half-space.  We  cmscuss  the  implications 
of  this  model  for  the  generation  of  S  waves  by  block  movement  in  the  vicinity 
of  an  underground  explosion.  \ 
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I.  SUMMARY 

The  research  performed  under  this  contract  during  the  period 
1  October  1978  through  31  March  1979  can  be  divided  into  three  main 
topics,  upper  mantle  studies,  propagatlonal  distortions  caused  by 
non-planar  structure  and  S-wave  generation  by  block  movements  near 
explosions . 

In  Section  II,  we  discuss  the  upper  mantle  structure  of  north¬ 
western  Eurasia  where  long-  and  short-period  WWSS  seismograms  from 
nuclear  explosions  in  the  Soviet  Union  have  been  incorporated  with 
published  apparent  velocity  data  to  derive  a  compressional  wave 
velocity  model  for  the  upper  mantle  beneath  northwest  Eurasia. 

The  waveforms  were  modeled  using  generalized  ray  theory  and  a  starting 
model  was  perturbed  by  trial  and  error  until  a  fit  was  achieved.  The 
seismograms  from  9  to  13“  exhibit  impulsive  first  arrivals  (P^) , 
implying  a  smooth,  positive  velocity  gradient  between  depths  of  60 
and  130  km.  A  consistent  pulse  arriving  about  2  seconds  after 
at  the  ranges  11  to  17*  is  modeled  by  positioning  a  low  velocity 
zone  between  150  and  200  km.  The  P^  phase  becomes  diffracted  beyond 
15°  although  there  is  considerable  evidence  of  lateral  variation  at  where 
this  geometric  cut-off  occurs.  The  model  is  relatively  smooth  from 
a  depth  of  200  km  down  to  420  km  where  a  5Z  Jump  in  velocity  produces 
a  triplication  from  15  to  23°.  The  observations  from  21  to  26°  clearly 
show  another  discontinuity  at  a  depth  of  675  km  with  a  4Z  change  in 
velocity . 

In  Section  III,  glorified  optics  has  been  used  to  compare  the 
waveform  complications  caused  by  two  and  three  dimensional  structures. 
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A  scheme  Is  developed  to  synthesize  seismic  Images,  in  which  the 
direct  rays  give  a  preliminary  feature  of  a  structure,  while  the  multiple 
bounces  provide  further  information  fir  inproving  and  refining  the 
n»del.  A  method  is  discussed  for  inverting  non-9tratif led  structures 
like  sedimentary  basins  and  subductlon  plates. 

In  Section  IV,  we  consider  the  problem  of  reflection  and 
refraction  of  purely  conpressional  waves  Incident  on  an  interface 
separating  Identical  solid  half-spaces  in  which  the  condition  of 
continuity  of  shear  displacement  at  the  boundary  is  generalized  to  one 
that  allows  slippage.  The  problem  is  solved  using  the  Cagniard-de  Hoop 
technique.  It  is  found  that  the  generation  of  reflected  P  and  S  waves, 
as  well  as  transmitted  S  waves,  Is  most  effective  In  the  case  of  a 
perfectly  unbounded  half-space.  We  discuss  the  Implications  of  this 
model  for  the  generation  of  S  waves  by  block  movement  in  the  vicinity 
of  an  underground  explosion. 


II.  UPPER  MANTLE  STRUCTURE  OF  NORTHWESTERN  EURASIA 
J.  W.  Given  and  D.  V.  Helmberger 

I.  INTRODUCTION 

There  is  considerable  interest  in  structure  of  the  upper  mantle  with 
respect  to  the  current  problems  in  geodynamics.  Although  there  seems  to 
be  broad  agreement  on  the  major  structural  features  of  the  mantle,  the 
various  models  proposed  vary  significantly  in  detail.  Important  questions 
vet  to  be  resolved  are:  (1)  how  deep  are  the  so-called  "400"  and  "600" 
km  discontinuities;  (2)  how  large  are  the  velocity  jumps;  (3)  how  deep 
in  the  mantle  do  velocity  variations  exist;  and  (4)  how  well  do  these 
variations  correlate  with  other  plate  properties. 

To  answer  these  difficult  seismological  questions  requires  more  model 
resolution  than  generally  provided  by  the  classical  methods,  namely  travel 
times  and  a  sparse  set  of  (dT/dA)  measurements,  see  Wiggins  [1969].  Some  of 
the  ambiguity  of  earth  models  determined  by  classical  means  can  be 
eliminated  by  the  use  of  short-period  synthetic  seismograms.  Comparing 
the  relative  amplitudes  of  observed  short-period  multiple  arrivals  to 
synthetics  allows  a  better  determination  of  triplication  points  as 
discussed  by  Helmberger  and  Wiggins  [1971].  Several  regionalized  models 
for  the  western  North  American  continent  have  been  obtained  using  this 
approach,  see  Wiggins  and  Helmberger  [1973]  and  Dey-Sarker  and  Wiggins  [1976] 
Unfortunately,  short-period  seismograms  suffer  from  waveform  instability 
in  that  neighboring  stations  in  many  situations  do  not  produce  similar 
pulses.  On  the  other  hand,  long-period  body  wave  studies  indicate  that 
waveshapes  recorded  on  the  long-period  World-Wide  Standard  Seismograph 
Network  (WWSSN)  are  remarkably  coherent,  see  Burdick  and  Helmberger  [1978] 
and  Helmberger  and  Burdick  [1979].  Thus,  an  acceptable  model  must  yield 
synthetics  which  agree  with  the  observed  travel  times,  the  dT/dA  measurements 
tin1  re  I  a i  I vc  timing  and  amplitudes  of  short-period  pulses,  and  the  long- 
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period  waveshapes.  Accurate  determinations  of  the  earth's  Green's 
functions,  or  the  broad-hand  step  responses  used  In  this  synthetic  modeling, 
for  various  regions  would  be  very  useful  In  quantitatively  answering  the 
questions  posed  above.  To  this  end,  a  data-set  compiled  from  Soviet 
explosl  ons  was  used  to  obtain  an  upper-mantle  model  which  will  be  used  as 
a  reference  model  In  comparing  other  regions. 

There  have  been  a  relatively  large  number  of  nuclear  explosions  In 
western  Russia  since  the  deployment  of  the  WWSSN.  Many  of  these  events 
were  well  recorded  In  northern  Europe  at  triplication  ranges,  see  Fig.  1. 

One  particularly  appealing  feature  of  this  data  set  Is  the  apparent 
absence  of  a  shadow  zone  near  9°  that  Is  common  in  regions  of  western 
United  States,  see  Romney  et  al.  [1962]  and  Helmberger  119  7 3b } .  A  typical 
seismogram  at  produced  bv  a  N.Z.  event  Is  displayed  In  Fig.  2.  Note  the 
sharp  F  and  S  phases  with  little  signs  of  anv  differential  attenuation. 

This  feature  Is  most  easily  explained  by  restricting  the  propagation  path 
to  be  entirely  In  the  lithosphere.  Thus,  we  have  the  prospect  of  determining 
a  complete  model  whereas  in  previous  studies,  the  serious  nonuniqueness  caused 
bv  the  low  velocity  zone,  (l.VZ)  made  the  absolute  depths  of  the  upper  mantle 
discontinuities  difficult  to  resolve. 

In  addition,  the  NORSAR  array  is  favorably  located  to  measure  the 
apparent  velocity  of  the  different  phases  in  the  P-wavetrain,  and,  thus,  measure 
(dT/diM.  Several  recent  investigators.  King  and  Galcagnite  [1976),  Masse 
and  Alexander  | 19741,  and  others,  have  used  these  measurements  to  study  the 
upper  mantle.  We  incorporated  their  apparent  velocity  information  into  our 
study  bv  adopting  KCA,  the  model  proposed  by  King  and  Calcagnile  [1976],  as 
our  starting  model.  The  initial  mode  1  was  modified  as  required  by  the  waveform 
data  but  the  dT/dA  curve  was  perturbed  as  little  as  possible. 


5 


II .  SYNTHETIC  SE IS  HOC RAMS  FOR  EXPLOSIVE  SOURCES 

The  computational  procedure  used  in  generating  synthetics  is  well 
established,  namely  one  assumes  that  a  synthetic  can  be  represented  by 
a  series  of  linear  operators: 

v(t)  “  s(t)*m  (t)*  a(t)*  r(t)*  i(t). 

s(t)  represents  the  source;  m(t),  the  mantle  response;  a(t)  an  attenuation 
operator;  r(t)  isthe  receiver  structure;  and  i(t)  is  the  instrument. 

We  calculate  the  source  time  function  by  the  procedure  discussed  in 
Burdick  and  Helmberger  [1979].  The  time  function  parameterization  is 
taken  from  Von  Saggarn  and  Blaodford  [1972} .  The  far-field  time  function  is 
written 

s(t)  -  *  ~  U-e"kt  (l+kt-B(kt)2)  H(t)*  C(t) 

o  dt 

where  ^  is  an  amplitude  factor,  C(t)  is  an  operator  representing  the 
effects  of  the  free  surface,  and  B  and  k  are  adjustable  parameters  controlling 
the  shape  of  the  time  function. 

The  parameter  B  controls  the  tradeoff  between  an  impulsive  and  step¬ 
like  time  function.  A  larger  value  of  B  means  that  the  source  is  more 
step-like,  a  smaller  value  of  B  indicates  a  more  impulsive  time  function. 

In  larger  nuclear  explosions  (i.e.,  for  modeling  long-period  seismograms) 
we  found  that  B  *  6  was  an  appropriate  value  to  adequately  reproduce  the 
teleseismlc  observations.  For  smaller  explosions  and  short-period 
seismograms  a  value  of  B  »  2  was  found  to  be  suitable. 

The  parameter  k  is  used  to  adjust  pulse  width  and  rise  time.  Larger 
explosions  are  modeled  by  assuming  k  ■  1,  small  explosions  are  modeled  by 
taking  k  »  S. 
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The  parameters  B  and  k  trade  off  with  the  operator  C(t)  which  models 
the  effects  of  the  free  surface.  The  shallow  source  depth  causes  the 
negative  reflection  from  the  free  surface  to  filter  out  the  low  frequencies 
in  the  t imef unct ion  mimicking  the  effect  of  both  B  and  k.  We  found  that 
P-pP  delays  of  0.9  sec  for  large  events  and  0.6  sec  for  small  events  were 
appropriate  for  our  purposes.  The  reflection  coefficient  at  the  free  surface 
was  taken  to  be  -0.8.  As  we  calculate  synthetics  for  different  ranges,  the 
trade  off  angle  of  P  and  pP  change.  As  a  result  the  pP-P  delay  and  the 
free  surface  reflection  coefficient  vary  with  the  distance  to  the  receiver. 

We  have  found  that  the  change  in  take-off  angle  has  a  small  effect  on  the 
waveforms,  not  significant  in  view  of  the  other  uncertainties  in  the  source 
operator.  The  operator  C(t)  also  includes  the  crustal  response  at  the  source. 
Anomalously  large  S  waves  (as  in  Fig.  2),  observable  out  to  15°  suggest 
that  there  may  be  some  crustal  contamination  in  sources  at  Novaya  Zemlya. 

We  could  find  no  evidence  in  the  long-period  P  waves  for  any  large  distortion 
due  to  source  structure  so  the  effect  has  been  ignored.  In  short-period 
seismograms  source  structure  is  Important  and  causes  a  significant  amount 
of  noise  which  is  difficult  to  interpret.  However  an  uncertainty  has  been 
introduced  which  should  be  considered  in  the  interpretation. 

At  ranges  where  the  arrivals  from  the  upper  mantle  are  separated 
in  time,  the  uncertainties  due  to  the  source  function  are  minimal.  At  ranges 
where  the  arrivals  from  the  upper  mantle  interfere,  the  errors  in  the 
source  function  can  strongly  affect  the  synthetics.  Of  particular  difficulty 
are  secondary  p-wave  arrivals  which  interfere  with  the  second  upswing  Of 
the  first  arrival.  We  tested  the  different  sources  to  check  the  validity  of 
our  conclusions  in  these  situations. 

The  effect  of  anelasticitv  are  approximated  by  a  Futterman  attenuation 
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operator  [Futterman,  1962;  Carpenter,  1966]  with  T/Q  -  1  for  all  ranges. 
Kennett  [1975]  has  pointed  out  the  errors  introduced  by  this  approximation 
when  low  Q  zones  are  present  in  the  upper  mantle.  The  assunption  was 
reevaluated  by  Burdick  and  Helmberger  [1978]  and  they  concluded  that  the 
relative  amplitudes  (waveforms)  were  not  affected  significantly.  We 
do  not  have  a  data  set  which  can  be  used  to  constrain  Q  so  we  will  follow 
the  technique  used  by  Burdick  and  Helmberger  [1978]. 

The  receiver  operator  models  the  effect  of  locating  the  receiver 
on  the  free  surface  and  the  effect  of  crustal  structure  beneath  the 
receiver.  The  scale  of  the  inhomo gene i ties  in  the  crust  are  often  small 
enough  to  be  ignored  on  the  long-period  seismograms,  see  Burdick  and 
Langston  [1977].  However,  crustal  structure  is  probably  the  most  important 
factor  contributing  to  the  inconsistency  observed  in  short-period  body 
waves  as  can  be  surmised  by  examining  the  horizontal  components  of  motion 
of  complicated  P-waveforms. 

The  other  operators  have  been  fixed  so  that  we  may  investigate  the 
operator  m(t),  the  upper  mantle  response.  The  earth  is  assumed  to  be 
laterally  homogeneous,  and  the  radial  velocity  structure  is  approximated 
by  layers  of  constant  velocity.  An  earth  flattening  transformation 
approximates  the  spherical  earth  by  a  plane  layered  medium.  The  Green's 
function  is  approximated  by  summing  a  finite  number  of  generalized  rays, 
the  individual  ray  response  is  calculated  using  a  Cagniard-de  Hoop 
algorithm  [see  Helmberger,  1973a  for  details].  To  maximize  computations 
the  ray  sum  is  truncated  after  including  all  of  the  primaries  for  a  time 
period.  In  ranges  where  diffraction  is  important,  in  shadow  zones  and  off 
the  ends  of  triplications,  multiple  reflections  were  added  to  test  the 
convergence  of  ray  sum.  The  accuracy  and  limitation  of  this  technique  is 
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discussed  by  Burdick  sud  Circuit  [1979]  and  will  not  be  reviewed  here  in 
the  Interest  of  brevity. 

At  the  ranges  and  for  the  structures  we  will  be  considering,  the 
phase  PP  Is  an  important  consideration.  An  anomalous  later  phase  can  be 
Identified  at  several  stations  12  to  15  sec  after  the  first  arrival  which 
we  tentatively  ldentlfi*  1  as  PP.  Preliminary  synthetic  models  of  PP 
Indicate  that  this  phase  could  be  large  depending  on  the  internal 
reflections  at  the  surface  near  the  halfway  point. 

III.  DETAILED  MODEL  DETERMINATION 

We  will  begin  this  section  by  presenting  our  results,  namely  model  K8, 
followed  bv  a  comparison  of  (dT/dA)  measurements  and  travel  times  with  the 
model  predictions.  Next,  we  will  give  a  detailed  comparison  of  the 
waveform  data  with  synthetics  generated  from  K8. 

The  starting  model  KCA  and  model  K8  are  shown  in  Fig.  3  and  their 
corresponding  (dT/dA)  values  displayed  in  Fig.  4.  There  are  two  important 
differences  between  the  models.  One  is  the  low  velocity  zone  at  150  km 
depth  in  K8.  The  other  is  the  steep  positive  velocity  gradient  between 
300  and  400  km  in  K8  where  KCA  has  a  very  small  gradient. 

Travel  time  curves  for  the  two  models  are  compared  in  Fig.  5. 

Plotted  with  the  calculated  curves  are  readings  from  the  ISC  for  the  nuclear 
explosions  in  Novaya  Zemlya  during  the  period  1966-1974. 

However,  the  inconsistency  of  observed  PP  from  station  to  station  makes 
this  phase  difficult  to  use  as  a  model  discriminant.  Upon  examination 
of  this  figure  it  should  be  noted  that  KCA  predicts  the  travel  times  at 
ranges  18  to  22®  and  23  to  25°  better  than  K8  although  diffracted  arrivals 
tend  to  prolong  the  A  branch  out  to  at  least  18°.  Furthermore  there  appears 


to  be  scattered  energy  arriving  between  the  Ad  and  Bd  arrivals 
which  can  be  attributed  to  high  frequency  scattering  phenomena.  It  i9 

plausible  that  diffraction  and  scattering  ran  extend  the  Ad  branch  as 

far  out  as  20°  making  the  K8  model  more  consistent  with  the  data. 

The  deviation  of  K8  from  the  data  at  around  23°  1b  harder  to 

justify.  In  this  case,  we  fit  the  waveform  data  at  the  expense  of  the 

travel  times.  It  is  possible  that  the  discrepancy  is  indicating  that  our 
model  is  too  simple  but  the  data  set  is  not  extensive  enough  to  resolve  this 
question. 

Crustal  structure  (0  to  40  km) 

The  crustal  structure  of  model  K8  was  taken  from  KCA  which  assumed 
a  single  40  km  layer  with  velocity  of  6.4  km/sec.  This  structure  was 
derived  from  a  synthesis  of  published  data  by  Der  and  Landisman  [1972]. 

The  vertical  travel  time  is  consistent  with  the  three  layer  model  derived 
by  Masse  and  Alexander  [1974]  and  is  slightly  shorter  than  the  more 
recently  propesed  model  by  Vinnik  et  al.  [1978].  There  is  undoubtedly 
some  lateral  variation  off  the  edges  of  the  east  Europear  platform  that 
are  larger  than  the  differences  between  these  models  and  we  should  expect 
to  see  at  least  a  second  of  scatter  in  travel  times  at  triplication  distance: 
see  Vinnik  et  al.  [1978],  Since  the  details  of  the  crustal  layering  has 
little  effect  on  the  synthetic  P-waveshapes  beyond  9°,  we  felt  Justified 
in  choosing  the  simplest  model  available. 

Lid  and  Low  Velocity  Zone  (40  to  300  km) 

The  profiles  of  long-  and  short-period  observations  which  were  used 
to  interpret  the  upper  300  km  are  shown  in  Fig.  6.  The  long  period  data 
is  exclusively  from  Novaya  Zemlya  sources,  the  short-period  observations 
are  from  Novaya  Zemlya  and  other,  smaller  explosions  to  the  south.  The 
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closest  rt  auah  If  short -period  data  are  at  14.5*.  Because  of  the  normal 
differences  In  gain  settings,  we  were  unable  to  obtain  a tmul taneoua 
recordings  of  any  one  event.  However,  the  multiplicity  of  both  large 
and  small  explosions  at  Novava  ?.emlya  corroborated  our  observations  of 
the  Important  features  In  the  data. 

The  I’-waves  from  ‘l  to  11.4“  are  relatively  simple  pulses  suggesting 
that  the  propagation  path  is  smoothly  varying.  The  large,  later  phase  at 
KBS,  11  sec  after  the  first  arrival.  Is  most  easily  explained  by  PI’;  however. 
It  Is  anomalously  large  at  this  particular  station  and  it  does  not 
appear  In  most  of  the  other  observations. 

Beyond  l  1.  1“  the  waveforms  begin  to  show  more  complexity,  the  long- 
and  short-period  setsmograma  showing  consistent  features  which  are 
Interpreted  as  upper  mantle  structure.  The  most  obvious  feature  on  these 
records  Is  the  onset  of  the  CD  branch,  the  reflection  from  the  400  km 
discontinuity.  Beyond  15“  this  phase  dominates  both  the  long-  and  short- 
period  seismograms.  Another  consistent  observation  is  the  small  long- 
period  first  arrival  at  14.5®  which  gtows  weaker  with  distance,  becoming 
almost  unobservable  on  the  seismogram  In  Fig.  b.  A  second  phase  become i 
apparent  on  the  long-period  data  at  UME  2-3  aec  after  the  firat  arrival. 

This  second  phase  Is  also  consistently  verified  on  the  short-period 
seismograms  although  It  is  of  a  variable  nature.  The  onset  of  the  second 
arrival,  and  the  dramatic  decay  of  the  amplitude  of  the  first  arrival  is 
crucial  to  the  argument  for  the  existence  of  a  low  velocity  zone  near 
21)0  km. 

We  perturbed  the  Initial  model  KCA,  to  fit  the  long-period  records  at 
UME  and  NUK  beyond  15“  while  trying  to  preserve  the  simple  nature  of  the 
waveforms  at  the  nearer  ranges.  Our  best  effort  at  matching  the  long- 
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period  features  Is  displayed  In  Kin.  7  where  we  have  Included  the  synthetics 
tor  model  Ki'A  for  comparison.  Ah  mentioned  earlier  the  convolution  ot 
the  source  function  given  In  the  right-hand  corner  with  the  derivative  of 
the  step  responses  vteld  the  synthetics.  Note  that  hy  Including  the  LVZ 
we  were  able  t o  produce  the  double  arrival  with  the  proper  decay  of 
the  Initial  arrival,  P(Ad)  with  range. 

The  corresponding  short-period  synthetics  are  displayed  In  Klg.  8. 

The  relative  timing  of  the  phases  In  the  synthetics  are  consistent  with  the 
data  as  well  as  some  of  the  relative  amplitudes. 

Transition  Region  ^ '00  to  800  Vm) 

The  seismograms  used  In  modeling  the  transition  region  are  particularly 
Interesting  as  can  he  surmised  hv  noting  the  intersections  of  the 
triplication  plot  a  In  Klg.  S. 

As  In  the  previous  section,  we  relied  heavily  on  the  long-period 
waveform  data  In  determining  the  branch  positions  and  the  resulting  model. 

The  synthetics  and  representative  waveforms  for  these  ranges  are 
displayed  In  Klg.  ‘J.  The  observed  wsveforms  are  somewhat  variable  In 
quality.  Kor  Instance,  the  K.0N  (21.1*)  observation  Is  the  only  seismogram 
at  that  range  from  an  explosion  at  the  southern  Novava  Zernlva  site  and  Is  a 
weak  recording.  On  the  other  hand,  NON  (21. *•*)  Is  from  the  northern 
site  where  large  shots  were  abundant,  some  examples  are  given  in  Klg.  10 
where  events  1  and  2  are  the  same  events  pictured  in  Klg.  7.  Thus,  we 
considered  tills  waveform  as  particularly  worthwhile  to  modal.  The  same 
situation  occurs  for  the  COP  (24. b*)  observation.  In  general,  these 
waveforms  are  extremely  difficult  to  model  using  a  trial  and  error  Inversion 
technique  because  each  record  Is  essentially  an  Interference  phenomenon  where 
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small  changes  In  the  source,  see  Fig.  10,  or  model,  see  step  responses,  can 
drastically  alter  the  waveform. 

One  of  the  principle  differences  between  K8  and  KCA  Is  the  gradient 
above  the  400  km  discontinuity.  Part  of  the  difference  Is  due  to  the  low 
velocity  zone  which  changes  the  bottoming  depth  of  rays  from  above  420 
km.  To  predict  correctly  the  relative  arrival  times  we  increased  the 
gradients  between  300  and  4-0  km  and  reduced  the  size  of  the  400  km 
discontinuity.  With  these  changes  we  were  able  to  model  the  interference 
at  21.9°  and  the  cut-off  of  the  long-period  AB  branch.  At  21.1°  the 
relative  arrival  times  of  the  phases  predicted  by  K8  fit  better  than  those 
calculated  for  KCA.  However  there  is  some  mismatch  in  the  relative 
amplitudes  with  the  arrival  from  above  the  420  km  discontinuity  being 
too  large.  The  choice  of  velocity  gradient  below  270  km  depends  crucially 
on  how  we  modeled  the  low  velocity  zone.  Any  errors  in  determining  the 
model  or  any  lateral  heterogeneity  in  the  upper  250  km  will  seriously 
affect  our  inversion  for  the  structure  between  280  and  400  km.  The 
uncertainties  in  K8  are  probably  greatest  in  this  depth  range. 

The  evidence  for  a  discontinuity  at  675  km  is  clear  in  King  and 
Calcagnile's  array  data  and  in  our  long-  and  short-period  observations. 

The  reflection  is  first  evident  at  KON,  21.1°,  and  forms  a  clear  distinct  pulse 
at  23.3°  (COP)  where  the  E-F  phase  is  the  second  pulse.  At  25.1°  AKU,  the 
E-F  arrival  is  first  and  the  shoulder  on  the  waveform  is  interpreted  to 
be  the  C-D  phase  moving  back.  The  best  data  showing  the  position  of  D  on 
the  E-D  branch  is  displayed  in  Fig.  11  where  Lhe  observations  obtained  from 
the  Semipalatln8k  test  site  are  compard  with  synthetics.  Due  to 
fortuitous  instrument  setting  near  26°  we  were  able  to  obtain  short-  and 
long-period  responses  for  the  same  explosions.  The  waveforms  appear  to  be 
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of  excellent  quality  and  are  reproducible  for  many  events.  A  careful 
examination  of  the  short-period  responses  at  thi9  particular  station 
indicates  that  it  is  relatively  transparent  as  judged  from  azimuthal 
and  horizontal  motion  tests  [see  Helmberger  and  Wiggins,  1971  for 
a  more  detailed  description].  For  this  reason,  we  have  spent  much 
effort  in  modeling  these  short-  and  long-period  observations  a9  closely 
a9  possible.  The  main  consequence  of  these  efforts  is  the  increased 
velocity  gradient  between  670  and  750  km,  see  Fig.  2,  which  increases 
the  amplitude  of  the  first  arrival.  Such  a  feature  is  also  consistent 
with  the  long-period  observations  at  AKU  (25.1°).  At  still  larger 
ranges,  we  see  no  compelling  evidence  of  either  the  (ED)  or  (AB) 
branch  although  the  short-period  observations  are  somewhat  perplexing 
with  respect  to  the  latter  branch,  see  Figs.  12  and  13. 

If  we  restrict  our  attention  to  only  NZ  events  we  see  relatively 
simple  seismograms  beyond  27°  which  are  in  agreement  with  our  model. 

On  the  other  hand,  Kazakh  events  tend  to  show  a  major  secondary  arrival 
out  to  larger  ranges  (labeled  by  an  arrow  in  Figs.  12  and  13).  Since 
King  and  Calcagnite  [1976]  used  mostly  Kazakh  events  we  can  see  why 
their  AB  branch  extends  to  larger  ranges.  This  feature  in  the 
synthetics  for  their  model  is  the  result  of  the  low  gradient  above 
420  km.  Thus,  this  difference  in  data  suggests  lateral  variations 
in  velocity  as  deep  as  400  km  with  respect  to  northern  and  southern 
Eurasia.  The  difficulties  with  this  interpretation  are  that  events 
from  Semipalatinsk  do  not  show  this  feature,  see  Fig.  11,  and  that 
many  Kazakh  events  do  not  show  this  branch,  for  example  see  KEV  (10) 
and  STU  (16).  It  should  also  be  noted  that  the  onset  of  the  AB  branch 
becomes  less  impulsive  with  distance  and  looks  more  like  the  onset  of 
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the  AB  branch  becomes  less  Impulsive  with  distance  and  looks  more  like 
the  onset  of  a  wavetraln,  see  Figs.  2k-2m  of  King  and  Calcagnlte  (1976]. 

Perhaps,  this  phenomenon  could  be  caused  by  multlpathlng  due  to 
small  scale  Irregularities.  It  would  be  particularly  interesting  to 
obtain  some  long-period  seismograms  that  sample  this  region  so  that  a 
comparison  similar  to  Fig.  14  could  be  made. 

IV.  DISCUSSION 

It  Is  interesting  to  compare  model  K8  with  models  proposed  for 
western  United  States  (WUS),  see  Fig.  15,  where  the  same  basic  modeling 
techniques  were  applied.  Models  HWA  and  HWB  were  derived  from  the 
(dT/dA)  measurements  (Johnson,  1967],  the  travel  times  from  the  NTS 
explosions  and  short-period  waveform  data  as  discussed  earlier.  The 
travel  times  and  (fit)  separations  between  the  triplication  arrivals 
for  HWB  and  HWA  are  substantially  different  between  13  and  19°,  where 
the  first  arrivals  appropriate  for  the  HWA  region  is  from  3  to  6  sec 
later  than  for  HWB.  Model  T/  was  constructed  to  fit  the  HWA  data  set 
as  well  as  the  long-period  waveforms  obtained  from  a  series  of  well- 
studied  earthquakes.  On  the  other  hand,  the  travel  times  and  (fit)'s 
for  the  model  HWB  are  compatible  with  those  predicted  by  K8.  Furthermore, 
a  comparison  of  the  short-period  seismograms  displayed  in  Fig.  6  are 
quite  similar  to  the  short-period  profile  presented  by  Helmberger  and 
Wiggins  (1971],  Fig.  2.  At  ranges  beyond  20°,  the  distinction 
between  the  regionalization  HWA  and  HWB  disappears  and  there  appears 
to  he  little  evidence  of  lateral  velocity  variations  below  300  km. 

Note  that  T7  and  K8  are  remarkably  similar  in  structure  below  this 

H 

depth.  Tiie  slight  off-sets  in  discontinuities  reflect  the  lid  and  LVZ 
disparity  which  is  worth  a  brief  review. 


15 


The  upper  200  km  of  T7  was  constructed  to  fit  the  average  travel 
times  and  amplitude  properties  of  WUS  where  there  appears  to  be  an 
effective  shadow  zone  along  some  profiles,  see  Helmberger  [1973b].  This 
effect  can  be  caused  by  a  LVZ  as  In  T7  but  It  could  be  caused  by  lateral 
changes  In  velocity  when  crossing  the  various  structural  provinces, 
for  instance  see  York  and  Helmberger  [1973].  Thus,  the  basic  vertical 
structure  of  HWA  or  K8  may  be  appropriate  for  some  portions  of 
northern  WUS  and  the  known  short-period  amplitude  pattern  caused  by 
horizon  change  in  structure.  Small  changes  in  the  velocity  structure 
of  the  Lid  and  LVZ  can  make  a  dramatic  change  in  synthetic  seismograms 
at  the  near-in  distances  of  8  to  19°  as  can  be  seen  by  the  following 
case  study. 

In  collecting  a  data  set  of  seismograms  one  invariably  finds 
an  odd  observation  which  appears  incompatible  in  travel  time  and  wave¬ 
shape  with  the  others,  such  is  the  case  with  the  NUR  record  displayed 
in  Fig.  16.  NUR  is  a  recording  of  a  southern  Novaya  Zemlya  event 
whereas  UME  was  produced  by  a  northern  Novaya  Zemlya  event,  see  Fig.  1. 

As  with  nearly  all  of  the  long-period  observations,  the  NUR  record  is 
duplicated  for  two  different  events  so  that  an  unusual  source  function 
does  not  appear  to  be  a  likely  explanation  for  the  difference  in  waveforms. 
Explaining  this  record  in  terms  of  a  perturbation  of  model  K8  proved 
to  be  quite  easy  and  enlightening.  By  slowing  the  first  arrival  by 
about  1  sec  relative  to  the  reflection  from  below  the  low  velocity  zone 
and  the  420  km  discontinuity,  we  can  reproduce  the  relative  timing  of 
the  seismograms.  The  result  is  model  K8'  shown  in  Fig.  16.  Note 
that  we  have  also  slightly  decreased  the  depth  to  the  low  velocity  zone. 
Another  characteristic  of  the  model  is  that  a  critically  reflected 


arrival  now  cones  from  below  the  low  velocity  zone;  i.e.,  the  shadow 
zone  ends  at  15*  rather  than  19*.  This  feature  substantially  Increases 
the  amplitude  of  the  second  arrival  relative  to  the  C-D  reflection. 

We  conclude  that  the  upper  mantle  in  northwestern  Eurasia  is  some¬ 
what  variable.  Model  K8  may  average  over  variations  as  large  as 
those  between  K8  and  K.88  so  that  the  detailed  features  of  K8  may  be 
considered  to  be  exemplary.  The  trade-offs  involved  in  modeling  low 
velocity  zones  are  well  known,  see  Dowling  and  Nuttli  [1964],  Helmberger 
[1973b].  To  more  fully  constrain  the  features  of  a  low  velocity  zone 
we  need  to  incorporate  longer  period  data  from  earthquake  sources. 

Fig.  18  complements  Figs.  !6  and  17  by  comparing  long-period  and  short- 
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period  synthetics  for  K8  and  K88.  The  long-period  seismograms  are  from 
an  earthquake  source  appropriate  for  the  March  23,  1978  Bermuda  event 
(Gordon  Stewart,  work  in  progress).  The  long-period  seismograms  are 
verv  similar,  hence  their  usefulness  in  constraining  the  overall 
features  of  the  model  such  as  the  average  gradients  and  the  sizes  of 
discontinuities.  The  short-period  data  are  more  sensitive  to  small  scale 
details.  The  need  to  consider  broad-band  inf ormat ion ,  long-  and  short- 
period  waveshapes  as  well  as  dT/dA  and  travel  time  data,  is  obvious  from 
this  series  of  figures. 

The  major  new  feature  in  model  K8  is  the  low  velocity  zone  and  it 
is  based  on  the  long-  and  short-period  observations  between  9  and  18°. 

It  is  not  possible  to  completely  rule  out  a  simple  discontinuity  at  220 
km  depth  as  has  been  suggested  by  Lehmann  [1964,  1967]  for  the  U.S. 
Several  models  including  such  a  feature  were  tested  but  none  were  found 
that  fit  the  data  as  well  as  K8.  Such  a  model  would  be  required  to 
have  a  nearly  constant  velocity  over  substantial  region  of  the  upper  40Q 
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km.  In  addition,  a  different  Q  .structure  than  implied  in  this  study 
would  be  necessary  to  predict  the  correct  amplitude  behavior  of 
the  travel  t i  me  branches.  As  our  knowledge  of  the  velocity  and 
Q  structure  of  the  earth  Increases  the  data  is,  of  course,  subject 
to  reinterpretation  under  other  working  hypotheses. 

It  is  mi  fort  unate  that  we  do  not  have  a  very  dense  coverage 
of  long-period  receivers  at  distances  of  19-2!)°.  The  high  frequency 
explosion  source  would  have  been  Ideal  for  elucidating  any  additional 
structure  between  the  two  m.i)or  discontinuities.  Models  HWB,  SMAK 
[Simpson  et  al.,  197bl,  and  ARC-TR  [Fukao,  1977]  all  have  an 
Inflection  in  the  velocity  structure  at  about  550  km.  Burdick 
and  Helmberger  [1978]  have  pointed  out  the  difficulties  in  observing 
this  feature:  the  short-period  data  have  too  much  resolving  power, 
the  long-period  earthquake  data  have  too  little.  We  found  no 
direct  evidence  to  Justify  Including  such  a  feature;  however,  considering 
the  discrepancies  in  the  observed  and  predicted  waveforms,  the  structure 
between  t  tie  two  dlscont inut le  is  remains  somewhat  uncertain. 

In  summary,  we  have  Incorporated  the  long-period  and  short-period 
waveforms  data  with  travel  time  and  apparent  velocity  measurements  to 
derive  the  model  K8  (Table  2)  for  Northwest  Eurasia.  The  significant 
features  of  the  model  are: 

(1)  A  low  velocity  /.one  between  depths  of  150  and  200  km; 

(2)  A  4.5"  velocity  Increase  at  420  km  depth; 

(3)  A  large  velocity  gradient  between  420  and  675  km; 

(4)  A  42  velocity  increase  at  675  km  depth. 
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Model  K8  is  similar  to  model  T7  proposed  recently  for  the  Western 
United  States;  however,  without  a  detailed  determination  of  the 
uppermost  velocity  structure,  any  discussion  of  lateral  heterogeneity 
is  premature. 
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Table  2.  Velocity  Model  K8.  Depths  Are  to  Top  of  Layer 


Velocity  Depth  Velocity  Depth  Velocity  Depth  Velocity  Depth 


6.39 

0 

8.392 

250 

9.572 

496 

10.955 

749 

6. AO 

1 

8.419 

260 

9.583 

506 

10.998 

759 

8.170 

40 

8.445 

270 

9.624 

517 

11.040 

770 

8.187 

50 

8.478 

280 

9.667 

528 

11.061 

783 

8.205 

60 

8.511 

290 

9.708 

539 

11.083 

796 

8.222 

70 

8.545 

300 

9.750 

550 

11.105 

808 

8.239 

80 

8.578 

310 

9.792 

560 

11.127 

821 

8.256 

90 

8.612 

320 

9.834 

570 

11.148 

834 

8.27A 

100 

8.645 

330 

9.876 

580 

11.170 

847 

8.291 

110 

8.678 

340 

9.818 

590 

11.192 

859 

8.308 

120 

8.712 

350 

9.960 

604 

11.213 

872 

8.325 

130 

8.745 

360 

10.003 

610 

11.235 

885 

8.3A3 

140 

8.764 

370 

10.045 

620 

11.256 

898 

8.250 

150 

8.783 

380 

10.090 

631 

11.278 

911 

8.180 

155 

8.803 

391 

10.135 

641 

11.300 

923 

8.050 

160 

8.822 

400 

10.180 

653 

11.321 

936 

8.0A0 

170 

8.841 

410 

10.225 

664 

11.343 

949 

8.150 

180 

9.24 

420 

10.660 

675 

11.365 

962 

8.270 

190 

9.292 

431 

10.702 

686 

11.385 

974 

8.287 

200 

9.333 

442 

10.744 

696 

11.408 

987 

8.305 

210 

9.375 

452 

10.789 

707 

8.323 

220 

9.415 

467 

10.829 

717 

8.340 

230 

9.458 

474 

10.891 

728 

8.366 

240 

9.50 

485 

10.913 

738 

i 


■  .V  .. 


23 


Table  1.  Nuclear  Explosions 


Event  No. 

Date 

Hr :Mln:S«c 

Lat  E 

Lon  N 

1 

08/28/72 

5:59.57 

73.30 

55.1 

Novaya  Zemlya  N 

2 

09/27/71 

5:59.55 

73.40 

55.1 

Novaya  Zemlya  N 

3 

11/02/74 

4:59.57 

70.80 

53.91 

Novaya  Zemlya  S 

4 

10/27/73 

6:59.57 

70.80 

54.20 

Novaya  Zemlya  S 

5 

09/12/73 

6:59.54 

73.30 

55.20 

Novaya  Zemlya  N 

6 

10/27/66 

5:58.00 

73.38 

54.62 

Novaya  Zemlya  S 

7 

07/10/71 

17:00.00 

64.2 

54.77 

8 

11/08/68 

10:02.05 

73.40 

54.90 

Novaya  Zemlya  N 

9 

08/14/74 

14:59.58 

68.91 

75.90 

10 

09/26/69 

6:59.56 

45.89 

42.47 

Kazakh  E 

11 

08/20/72 

2:59.58 

49.46 

48.18 

Kazakh  N 

12 

10/22/71 

6:02.57 

51.51 

54.54 

13 

09/27/73 

6:59.58 

70.76 

53.87 

Novaya  Zemlya  S 

14 

07/01/68 

4:02.02 

47.92 

47.95 

Kazakh  N 

15 

10/03/72 

8:59.58 

46.85 

45.01 

Kazakh  E 

16 

12/22/71 

6:39.56 

47.87 

48.22 

Kazakh  N 

17 

12/06/69 

7:02.59 

43.83 

54.78 

Kazakh  W 

18 

12/12/71 

7:00.57 

43.85 

54.77 

Kazakh  W 

19 

12/23/70 

7:00  57 

43.83 

54.85 

Kazakh  W 

20 

02/13/66 

4:51.58 

49.82 

78.13 

Semlpalatlnsk 

21 

07/23/73 

1  -.22.58 

49.99 

78.85 

Semipalatinsk 

22 

12/18/66 

4:57.58 

49.93 

77.73 

Semipalat lnsk 

23 

02/26/67 

3:57.58 

49.78 

78.12 

Semlpalatlnsk 
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Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


Ft  gvire 

Figure 

Fi gure 


Capt Iona 

1.  Map  Indicating  the  location*  of  aource  events  and  WWSSN  seismographs 
used  in  this  study.  Information  about  the  event  numbers  is  given 

in  Table  1. 

2.  Umg-period,  vertical  component  seismogram,  KEV  from  Novava  Zemlva 
(8/28/72)  recorded  at  A  ■  9.4°.  The  s-phase  is  much  larger  than 
the  p-phase  and  has  virtually  the  same  frequency  content. 

3.  Model  KCA  and  K8. 

4.  Comparison  of  the  (dT/dA)  measurements  reported  by  King  and 
Calcagnlte  [1976]  and  K8  predictions. 

5.  Comparison  of  predicted  travel  time  curve  for  K8  and  KCA  and 
travel  t imes  of  Novaya  Zemlva  explosions  as  reported  by  the  ISC. 

The  letters  on  the  K8  travel  time  plot  label  the  branches  for 
further  reference  in  the  text.  Solid  lines  Indicate  direct 

ray  theoretical  arrivals,  dashed  lines  indicate  the  prolongation 
of  certain  branches  by  diffraction  due  to  a  low  velocity  zone. 

6.  Representative  seismograms  obtained  from  Novaya  Zemlya  explosions 
displaying  various  arrivals.  The  lines  are  keyed  to  the  triplication 
plots  given  in  Fig.  4.  The  number  following  the  station  label 
indicates  the  event  (Table  1)  where  1  and  2  are  from  the  northern  site. 

7.  Comparison  of  synthetics  with  the  long-period  waveform  data. 

8.  Short-period  synthetics  for  KCA  and  K8  which  can  be  compared  with 
the  data  In  Fig.  6. 


i  _ _ ... 


t _ 


Figure  9.  Comparison  of  the  observed  waveforms  with  synthetics.  The  K8 

synthetics  were  produced  by  a  convolution  of  the  source  function 

displayed  in  the  right-hand  corner  with  the  derivative  of  the 
step  responses. 

Figure  10.  Comparison  of  synthetic  (K8)  with  observations  from  various  events 
at  the  Northern  Novaya  Zemlya  test  site  indicating  source 
variabi lity . 

Figure  11.  Observations  at  Tabriz,  Iran  from  nuclear  explslons  in  Semipalat insk . 

The  data  clearly  show  evidence  for  the  675  km  discontinuity  and  the 
position  of  the  end  of  the  C-D  branch. 

Figure  12.  Theoretical  seismograms  from  K8  fit  the  short-period  seismograms  as 
well  as  KCA.  The  inconsistency  of  the  records  between  24  and  25° 
should  be  noted.  The  large  later  phase  indicated  by  the  arrows  is 
similar  to  the  arrival  of  King  and  Calcagnlte  [19761  interpreted 
as  the  reflection  from  the  420  km  discontinuity. 

Figure  13.  The  figure  compares  theoretical  and  observed  seismograms  from  25 
to  29°.  Again  note  the  occasional  appearance  of  the  reflection 
from  the  420  km  discontinuity.  This  appears  to  be  a  regional 
phenomena,  observed  only  on  sources  in  Kazakh  and  receivers  in 
Northwest  Europe. 

Figure  14.  The  figure  shows  the  theoretical  and  observed  seismograms  of  Novaya 
Zemlya  explosions  recorded  at  STU.  Model  KCA  still  predicts  a  large 
reflection  from  the  420  km  discontinuity  which  is  not  evident  in 


the  data. 
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Figure  15.  Comparison  of  model  K8  with  models  HWA,  HWB,  and  T7  derived  from 
WUS  observations. 

Figure  16.  Seismograms  from  Northern  Novaya  Zemlya  to  UME  and  from  Southern 
Novava  Zemlya  to  NUR.  These  seismograms  are  along  slightly 
different  azimuths.  The  difference  can  be  explained  by  variation 
In  the  velocity  structure  above  180  km  (see  Fig.  17). 

Figure  17.  Model  K8*  which  was  derived  from  K8  to  explain  the  variation  at 
15. 6°,  see  Fig.  16. 

Figure  lb.  Short-  and  long-period  synthetics  for  K8  and  K8'.  The  long-period 
source  function  Is  from  the  March  23,  1978  Bermuda  earthquake 
(G.  Stewart,  personal  communication). 
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III.  TELESEISMIC  PROSPECTING  OF  NON-STRATI  FI  ED  RECEIVER  STRUCTURES 
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By  Tai-Lln  Hong 

I.  INTRODUCTION 

Seismic  prospecting  is  usually  conducted  in  two  ways.  The  first 
is  the  method  of  zero-offset  seismic  sections.  It  is  good  for  fine 
anomalous  structures  like  the  petroleum  reservoir  but  is  subject  to 
shallow  penetration  and  a  narrow  range  of  scope.  The  second  is  the 
method  of  refraction  profiles.  It  is  good  for  structures  whose  inter¬ 
faces  are  planar  but  not  able  to  find  curved  anomalies. 

In  this  paper  another  approach  using  teleseismic  body  waves  as 
sources  is  discussed.  This  approach  has  two  merits.  First,  it  car. 
sense  the  structures  in  the  depth  since  the  wave  comes  steeply.  Second¬ 
ly,  since  the  incident  wavefront  is  planar,  the  resultant  wave  forms  in 
the  receiver  structure  have  the  minimum  feature  of  dispersion;  therefore, 
the  high-frequency  approximations  can  be  used  to  recover  the  informa¬ 
tion  contained  in  the  wave  forms. 

As  a  matter  of  fact,  this  approach  has  been  discussed  for  the 
cases  of  planar  structures  by  Burdick  and  Langston  (1977)  and  Langston 
(1977)  among  others.  Therefore,  this  paper  will  stress  on  the  structures 


with  curved  boundaries. 
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II.  THE  PATH  OF  MULTIPLE  BOUNCES  IN  THREE-DIMENSIONAL  STRUCTURES 

In  a  three-dimensional  structure,  a  ray  can  be  bent  to  a  quite 
different  azimuthal  direction  after  several  bounces.  This  effect  is 
sensitive  to  the  structure  geometry  as  well  as  the  incident  direction 
and  results  In  various  waveform  distortions. 

In  this  section  we  will  illustrate  this  effect  through  several 
cases  in  which  the  structures  are  of  the  type  of  a  sedimentary  basin 
over  a  half-space.  The  S-velocitv,  P-velocity,  and  density  of  the 
material  In  the  basin  are  assumed  to  be  0.7  km/sec,  1  km/sec,  and 
2  g/cm*.  Those  In  the  half-space  are  assumed  to  be  3.5  km/sec,  5  km/sec 
and  3.3  g/em*. 

The  geometry  of  the  interface  can  be  expressed  as 

*  -  d  +  |  [i  -  Cos  2»(-x--+-  ^  ---w/2>  ]  (Eq.  1) 

where  x,  v,  and  z  are  the  coordinates  in  a  master  frame,  as  shown  in 
in  Fig.  1-a. 

c  and  d  are  the  parameters  related  to  the  thickness  of  the  basin 
w  is  the  parameter  related  to  the  width  of  the  basin;  it  is  a 
function  of  x  and  v  as  follows 

x‘  +  y* 

y  a  y  — — - d. - # 

°  (1-A2)  x2  +  y2 

The  incident  direction  of  the  plane  wave  is  expressed  in  terms  of 

the  azimuth  $  and  latitude  0  in  the  master  frame  as  follows: 
o  o 

0  =  (v  ,  v  ,  v  )  =  (sinB  eos4>  ,  sin0  simp  ,  cos0  ). 
x  y  z  o  o  o  o’  o 

For  the  convenience  of  the  scanning  procedure  which  will  be  illus¬ 
trated  In  the  next  section,  we  need  a  scanning  frame  x*~  y'-z'.  A 
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reference  point  (xq,  yQ,  zq)  Is  set  somewhere  below  the  receiver  struc¬ 
ture.  The  scanning  frame  Is  formed  by  translating  the  master  frame  to 
(xq,  yQ,  z^)  and  rotating  it  by  the  azimuthal  angle  $o  and  the  latitude 
angle  0q.  Therefore,  the  z'-axis  is  parallel  to  the  incident  direction 
of  the  plane  wave.  Then,  we  pick  an  azimuthal  direction  in  the  scanning 
frame  to  set  up  a  baseline  along  which  the  scanning  rays  are  sent  from 
the  plane  z'“  0. 

In  each  of  the  following  cases,  we  will  send  a  group  of  rays  of 
S  phase  along  the  9ame  baseline.  These  rays  are  transmitted  into  the 
basin,  reflected  twice  on  the  interface,  and  terminated  on  the  surface. 

Case  1 

The  parameters  for  the  geometry  of  the  interface  are: 

d  *  -1  km,  c  =  -5  km,  w  =  50  km,  and  A  =  0. 

The  incident  direction  of  the  plane  wave  is: 

<J>  =0.2  radians,  9  “  0.2  radians, 

o  o 

The  scanning  reference  point  (x^,  yQ ,  zq)  is  set  at  (0,  0,  -20  km). 

The  scanning  baseline  is  set  at  the  azimuth  0°  in  the  scanning  frame. 

Notice  that  the  structure  here  is  cyl indrically  symmetric.  Since 
this  group  of  rays  are  sent  on  a  plane  which  contains  the  axis  of 
symmetric,  there  is  no  azimuthal  scattering,  as  shown  in  Fig.  1-a. 

Case  2 

The  structure  and  the  incident  direction  of  the  plane  wave  are 
exactly  the  same  as  those  in  Case  1. 

The  scanning  reference  point  (xq,  yQ,  zq)  is  set  at  (-5  km,  -5  km,  -20  km). 

The  scanning  baseline  is  set  at  the  azimuth  0°  in  the  scanning  frame. 

Now,  this  group  of  rays  is  no  longer  on  a  plane  which  contains  the 
axis  of  symmetry;  therefore,  these  rays  are  scattered  azimuthally  as  shown 


in  Fig.  1-b.  We  fan  see  from  the  horizontal  projection  that  the 
sharpest  bending  of  the  scanning  rays  occurs  when  the  rays  are 
transmitted  into  the  basin,  because  of  the  large  velocity  contrast 
at  the  Interface. 


Case  1 

The  structure  is  exactly  the  same  as  that  in  Case  1. 

The  incident  direction  of  the  plane  wave  is 
$o  =  1.5708  radian,  0O  »  0.2  radian. 

In  other  words,  the  plane  wave  is  propagating  upward  along  the 
z-axis  and  forward  along  the  y-axis. 

The  scanning  reference  point  (xq,  y^ ,  zq)  is  set  at  (0,  0,  -20  km). 
The  scanning  baseline  is  set  at  the  azimuth  -1.5708  radian  in  the 
scanning  frame.  From  the  horizontal  projection,  we  can  see  that  the  rays, 
advancing  in  the  y-direction,  are  turning  around  after  the  transmission. 

At  the  edge  of  the  basin,  the  horizontal  projections  of  the  ray  segments 
in  the  basin  are  very  short  because  of  their  steep  paths. 

Case  A 

In  this  case,  everything  except  the  geometry  of  the  interface  is  the 
same  as  that  in  Case  1. 

The  parameters  for  the  geometry  of  the  interface  are: 
d  =  -1  km,  c  =  -5  km,  w^  =  50  km,  A  »  0.5. 

The  contours  of  this  interface  are  elliptical  as  shown  in  Fig.  1-a. 
Also,  notice  that  the  cross  section  of  this  basin  on  the  plane  y  *  0  is 
wider  than  those  in  the  previous  cases. 
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III.  THE  SCANN1N0  PROCEDURE 


Because  a  multiple  bounce  often  varies  its  azimuthal  direction 
along  its  path,  it  is  not  trivial  to  predict  the  starting  azimuth  for 
a  geometric  ray  leading  to  a  given  receiver  location.  A  scanning 
procedure  is  suggested.  The  trial  rays  are  sent  in  systematic  ways  to 
find  the  approximate  direction  of  a  geometric  path.  Then,  we  can 
take  a  fine  tuning  procedure  to  obtain  the  required  accuracy. 

In  this  section  we  will  Illustrate  the  scanning  procedure  for  a 
case  in  which  the  structure  is  of  the  type  of  a  two-dimensional 
sedimentary  basin  over  a  half-space.  The  materials  in  the  basin  and 
in  the  half-space  are  the  same  as  those  in  the  cases  of  Section  II. 

The  analytical  representation  of  the  geometry  of  the  interface  is  as 
follows : 


Cos 


2n(x  -  w/2)  J 


(Eq.  2) 


where  x,  z  are  the  cordinates  shown 

d  and  c  are  the  parameters  related  to  the  thickness  of  the  basin 
w  is  the  parameter  related  to  the  width  of  the  basin. 


In  this  case,  d,  c,  and  w  are  assumed  to  be  -1  km,  -5  km,  and 
50  km,  respectively. 

The  incident  direction  of  the  plane  wave  is: 

$o  ”  1.5708,  fl  =  0.31416  radians. 

Notice  that  this  plane  wave  is  propagating  not  only  upward  along 
the  z-axis  but  also  forward  along  the  y-axis.  T'nere fore,  although  the  structure 
geometry  is  not  varying  along  the  y-axis,  this  incident  direction  makes 
this  problem  actually  a  three-dimensional  problem. 

The  scanning  reference  point  (x  y  ,  z^)  is  set  at  (0,  0,  -20  km). 
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In  Fig.  2  are  shown  the  horizontal  projections  of  eight  groups  of 
scanning  ravs.  They  are  plotted  in  Fig.  2-a  through  2-h.  The  azimuth 
angle  of  the  baseline  of  each  group  in  the  scanning  frame  is  -2.0416, 
-1.7916,  -1.5416,  -1.0416,  -0.3,  -0.2,  -0.1,  0  radians,  respectively. 

The  scanning  rays  illustrated  here  are  all  in  S  mode.  They  are 
transmitted  into  the  basin,  reflected  four  times  on  the  interface,  and 
terminated  on  the  surface. 

In  Fig.  2-a  through  2-d,  the  patterns  of  azimuthal  scattering  are 
similar.  The  projections  of  the  ray  segments  in  the  half-space  are 
neatly  lined  up.  The  paths  of  the  first  ray  segments  in  the  basin  are 
verv  steep  because  of  the  large  velocity  contrast.  Thus  their  hori¬ 
zontal  projections  are  very  short.  The  ray  segments  in  the  layer  tend 
to  propagate  along  the  x-direction  since  the  interface  is  curved  along 
that  direction. 

In  Fig.  2-e  through  2-g,  we  can  see  more  clearly  how  these  ray 
segments  in  the  basin  are  bouncing  back  and  forth  near  the  cental  region 
of  the  basin. 

By  tracking  the  location  where  each  of  these  scanning  rays  ter¬ 
minates,  we  can  eventually  find  a  ray  tube  containing  the  receiver  for 
the  following  procedure  of  fine  tuning. 


IV.  WAVK  FORM  SYNTHKS I S 


The  glorified  opt  1  cs  ( see  Hong  and  Helmberger,  1978;  and  Hong, 
19/9)  Is  used  to  synthesize  tin*  wave  fonns.  It  Is  a  hint)  frequency 
approximation.  Ilir  reflection-transmission  effects  are  evaluated  In 
every  local  t rank*  formed  bv  the  normal  and  t augent 1 al  plane  to  a 
boundary.  The  spreading  effects  are  evaluated  globally  through 
Infinitesimal  rav  tubes.  The  displacement ,  u,  can  be  simply  written 
as  follows: 

u  -  V  7  fit  -  t(i)  Re  +  f*(t  -  I())  I  m  (,S'‘J)  j  (Eq.  1) 

n-1 

where  N  is  the  total  number  oi  revs  considered 

t  is  the  traveling  time  ol  each  ray  from  the  baseline  to  the 
rece i ve  r 

.Vis  the  spreading  effect 

'J  Includes  tlie  strength  ot  the  potential  tield  of  the  Incident 
plane  wave ,  t he  rei lect ion-t ransmission  effect,  and  the 
receiver  response 

f  is  the  convolution  ot  the  source  t  i  n>e  function  and  t  ho  instru¬ 
ment  response 

f*  is  the  Hilbert  transformation  of  I 

For  evaluat  ing.S'  ,  we  send  a  cylindrical  ray  tube  from  the  plane 
7.'  -  0  in  the  scanning  frame,  we  take  the  square  root  of  the  ratio  of 
cross  section  In  the  half-space  to  the  projected  cross  section  at  the 
receiver,  and  determine  the  9()1’  phase  shift  bv  inspecting  the  rotation 
of  the  element  rays  ot  the  tube.  The  rav  tube  Is  supposed  to  be 
Infinitesimal;  the  limiting  procedure  Is  taken  numerically. 


it 
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Kot  evaluat  iug'J  .  we  tlrst  Iltul  the  potential  vector  for  the 
incident  S  wave  oi  t  ho  scalar  potential  ^  for  tho  incident  I’  wave,  at,w 
trmisler  it  in  local  frames  to  obtain  tho  transmitted  or  reflected 
potentials.  rho  potential  ticld  reaching  tho  receiver  is  converted 
to  tho  dl  splnccment  field. 

A  more  detailed  description  ot  the  glorified  optics  for  three- 
dimensional  problems  can  ho  found  in  Hong's  (1979)  paper. 

The  nature  ot  wave-form  compl loat ion  in  some  two-dimensional 
structures  has  boon  discussed  lit  Hong  and  Helmberger's  paper  ll**7fO. 

Now,  wo  1  ike  to  see,  I  turn  a  case  ot  a  two  dimensional  structure  to 
a  ease  ot  a  t  lit  ee-d  linens  l  onal  si  meturo,  wli.it  kinds  of  further  com¬ 
plications  can  bo  added  to  tho  wave  torms.  Nilotic  speaking,  they 
can  bo  attributed  to  two  reasons.  First,  there  ate  more  rays  arriving 
from  various  a ••  i nmt ha  1  directions.  This  can  be  easily  conceived  hv 
looking  at  Fig.  i.  Secondly,  the  spreading  effect  is  sensitive  to  the 
curvature  in  azimuthal  directions.  To  illustrate  this,  we  consider  a 
two-dimensional  structure  like  that  in  the  case  ot  Section  111  and  a 
three-dimensional  structure  like  that  in  Case  1  ot  Section  II.  In  each 
case  we  examine  a  ruv  tube  surrounding  the  geometric  path  ot  a  multiple 
bounce  which  is  transmitted  into  the  basin,  reflected  three  times  on 
the  interface,  and  teimluated  at  t he  receiver  fin  km,  0,  0).  The  vert i 
cal  pro  lections  and  horizontal  protections  ot  these  ray  tubes  are 
almost  identical.  But  the  horizontal  project  ions  show  the  difference. 

We  notice  that  the  r.iv  tube  in  the  three-dimensional  structures  is 
twisted.  Such  i av  tubes  ot  large  size  are  lust  used  for  the  purpose  ot 
illustration.  Foi  computing  the  spreading  effect,  much  finet  tubes 


shoul d  be  used . 


q 


In  the  rem.it nt nit  part  of  this  section,  we  will  synthesize  and 
compare  the  displacement  wave  forms  for  three  groups  of  cases  to  show 
how  much  the  wave  forms  can  be  sensitive  to  the  azimuthal  variation  of 
the  geometry  of  a  structure.  Ricker's  wavelets  are  used  as  fit),  i.e.. 


fit) 


(Kq.  4) 


where 

«  .  (v  -  Jp) 2 

1  2 
(TP  -  ) 

In  each  of  the  following  cases,  we  use  two  flt)'s  with  Tp  equal 
to  18.3  sec  and  1  sec,  respectively.  Not  Ice  that  the  difference  between 
the  long  period  response  and  t lie  short  period  response  is  mainly  due 
to  the  traveling  time  of  each  arrival  instead  of  the  Intrinsic  frequen¬ 
cy  dependence.  The  latter  is  at  most  90°  phase  shift  In  the  high 
frequency  approx imat ion. 

1)  (.roup  A 

In  this  group  we  will  compare  the  responses  of  a  two-dimensional 
and  a  three-dimensional  structure  to  the  same  incident  S  plane  wave. 

The  structures  are  of  the  tvpe  of  a  sedimentary  basin  over  a  half-space. 
The  material  structure  is  the  same  as  that  described  in  Section  11.  The 
incident  direction  of  the  plane  wave  is  assumed  to  be  vertical.  The 
displacement  field  u"  of  the  incident  wave  is  assumed  to  be: 


u 


o 


10,  1 .  0) . 


Case  1 


The  structure  is  two-dimensional .  The  representation  form  for  the 


geometry  ol  the  interface  Is  (Kq.  2)  with  the  following  parameters: 
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J  *  -1  km,  c  •  -S  km,  w  «  SO  km. 

o 

Its  cross  section  on  the  plane  y  ■  0  and  Its  contours  are  plotted  In 
KIr.  4.  The  receiver  location  (x.  y,  z)  is  assumed  to  be  (0.,  0.,  0.,). 
The  scanning  reference  point  (x,  y^,  7.q)  is  set  at  (0,  0,  -20  km). 

The  traveling  timers',  and  ‘J  of  each  arrival  are  listed  in  Table  1. 

The  displacement  wave  forms  are  plotted  in  Fig.  S-a . 

Case  2 

Tlte  structure  is  three-dimensional.  The  analytical  representation 
of  the  geometry  of  the  interface  is  of  the  following  form. 


z  -  d 


II 


,,  2ir(x‘  + 
C  o  s 


w/2) 


(F.q.  5) 


where 


w 


w 

o 


fx'  +  v '  ) 
x-  +il  -A-  )v' 


The  parameters  are  assumed  as  follows: 

d  “-1  km,  c  -  -5  km,  -  SO  km,  A  -  O.S. 

Its  cross  section  on  the  plane  v  ■  0  and  its  contours  are  plotted  in 
Fig.  b.  The  receiver  and  the  scanning  reference  point  are  the  same  as 
those  in  the  last  ease. 

The  traveling  time, A,  and  U  of  each  arrival  are  listed  in  Table  2. 
By  comparing  the  ravs  No.  1  to  11  in  Table  1  and  those  in  Table  2,  we 
can  see  these  ravs  art’  taking  the  same  geometric  path  in  the  two  struc¬ 
tures,  but  their  spreading  effects  are  different.  We  can  see  the  higher 
concentrations  in  the  three-dimensional  structure.  Also,  we  notice  that 
in  the  t  h  reo-d  i  mens  i  on  a  1  structure  there  are  more  ravs  which  /ire  pro¬ 
pagating  on  the  plane  x  »  0,  Those  rays  are  listed  in  Table  2  as 


rays  No.  12  to  17.  However,  their  contribution  to  the  wave-form  distri¬ 
bution  is  slight.  The  reason  is  that  those  rays  are  reflected  as  SV 
waves  in  the  local  frames  and  thus  lose  more  energy  than  if  they  were 
local  SH  waves. 

The  displacement  wave  forms  of  this  case  are  plotted  in  Fig.  5-b . 
ii)  Group  B 

In  Case  2  of  the  last  group,  the  rays  No.  12  to  No.  17  do  not 
contribute  to  the  wave-form  distortion  significantly.  Now,  assume  the 
rays  are  in  P  mode;  those  rays  will  then  become  relatively  important. 

To  make  the  problem  simpler,  we  assume  the  materials  are  fluid.  The 

3 

P-velocity  and  density  are  assumed  to  be  0.7  km/sec  and  2  g/cm  in 

3 

the  basin  and  to  be  3.5  km/sec  and  3.3  g/cm  in  the  half-space.  The 

plane  wave  is  assumed  to  be  vertically  incident  with  the  displacement 

field  x  as  follows, 
o 

O  /  O  O  O  v  /  r\  r\  i\ 

'  =  ^Ux’  Uy ’  Uz)  =  (°»  °*  l) • 

Here  follow  Case  3  and  Case  4,  corresponding  to  Case  1  and  Case  2, 
respectively. 

Case  3 

The  structure  geometry,  the  receiver  location,  and  the  scanning  reference 
point  are  the  same  as  those  in  Case  1.  The  arriving  time,«S,  and 
are  listed  in  Table  3.  The  displacement  wave  forms  are  plotted  in 
Fig.  5-c. 

Case  4 

The  structure  geometry ,  the  receiver  location,  and  the  scanning  reference 
point  are  the  same  as  those  in  Case  2.  The  traveling  time.iS  ,  and  ^ 

.ire  listed  in  Table  4.  The  displacement  wave  forms  are  plotted  in  Fig.  5-d. 
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Now,  we  can  see  from  Fig.  5-d  that  rays  No.  12  and  No.  13,  which 
are  propagating  on  the  plane  x  =  0,  result  in  a  large  spike  which  does 
not  appear  in  Fig.  5-c. 

iii)  Group  C 

In  this  group  we  will  examine  the  response  of  a  two-dimensional 
structure  to  the  plane  waves  incident  in  different  azimuthal  directions. 
Hie  structure  used  in  this  group  is  the  same  as  that  in  Case  1  of  this 
section. 

Case  5 

The  incident  direction  of  the  plane  S  wave  is 
4>o  *  0,  0  =  0.31416  radians. 

Its  displacement  field  isu°=(u°,  u°,  u°)  =  (0,  1,  0). 

**  x  y  z 

The  receiver  location  (x,  y,  z)  is  set  at  (0,  0,  0). 

The  scanning  reference  point  (xq,  y^,  z  )  is  set  at  (0,  0,  -20  km). 

Notice  that  the  geometric  ray  paths  are  all  on  the  plane  y  =  0. 

The  rays  are  all  SH  waves  in  local  frames.  The  effects  of  the  con¬ 
centration  and  triplication  are  similar  to  those  discussed  in  the 
previous  cases  of  two-dimensional  structures.  The  traveling  time,iS, 
and  3  of  each  arrival  are  listed  in  Table  5.  The  displacement  wave 
forms  are  plotted  in  Fig.  5-e. 

Case  6 

The  incident  direction  of  the  plane  S  wave  is 

=  1.3708  radians,  0  = 

o  o 


0 . 31416  radians . 
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The  displacement  field  u  of  the  Incident  wave  is 

V  “  <ux-  Vly.  "z>  -  O.  0,  0). 

The  locations  of  the  receiver  and  the  scanning  reference  point 
are  the  same  as  those  in  the  last  case. 

Notice  that  the  geometric  ray  paths  are  all  on  the  plane  x  ■  0. 
The  rays  are  all  SH  waves  in  local  frames.  There  is  no  triplication 
because  the  cross  section  of  the  structure  on  the  plane  x  -  0  is  flat. 
However,  the  effects  of  concentration  3nd  phase  shift,  caused  bv  the 


curvature  ot  the  interface,  are  obviously  showing  up  in  Fig.  5-f.  The 
traveling  timc,.V,  and  'J  of  each  arrival  are  listed  in  Table  6. 
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V.  SEISMIC  IMAGE  OF  A  STRUCTURE  PRODUCED  BY  AN  INCIDENT  PLANE  WAVE 

There  are  many  parameters  contributing  to  the  wave-form  distortion 
in  a  receiver  structure.  It  is  very  difficult,  if  not  impossible,  to 
resolve  the  effect  caused  by  each  of  those  parameters  from  the  record 
at  a  single  station.  However,  if  we  are  running  a  seismic  array  in  a 
receiver  region,  there  could  be  enough  constraints  in  the  data  to 
uniquely  conclude  some  features  of  the  structure. 

For  this  purpose  we  developed  a  scheme  to  simultaneously  synthesize 
the  wave  form  at  each  station  of  an  array.  According  to  the  horizontal 
locations  of  those  stations,  the  seismograms  are  lined  up  to  simulate 
the  vertical  cross  sections  of  the  receiver  structures.  Such  an 
assembly  is  called  a  synthesized  seismic  Image  of  a  structure. 

Here  we  like  to  take  some  two-dimensional  structures  as  examples 
to  show  how  efficiently  the  seismic  Images  can  characterize  the  sub- 
surficial  structures.  These  structures  are  of  the  type  of  a  soft 
layer  over  a  half-space  with  a  curved  interface.  The  material  struc¬ 
ture  in  each  example  is  assumed  to  be  the  same  as  that  in  the  cases  of 
Section  II.  Also,  we  assume  the  Incident  direction  of  the  plane  wave 
in  each  example  is  consistent  with  the  two-dimensional  feature  of  the 
structure.  In  order  to  make  the  illustration  simpler,  we  will  consider 
only  the  direct  ray  and  the  first  two  multiple  bounces. 

i)  Example  1 

The  structure  Is  the  same  as  that  in  Case  1  of  Section  IV.  The 
direction  and  the  displacement  field  u^  of  the  incident  plane  S  wave 
are  also  the  same  as  those  in  Case  1  of  Section  IV. 

The  receivers  are  distributed  uniformly  on  the  surface  along  the 


x-axls  with  the  interval  of  1  km. 


1 
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The  Ricker's  wavelet  with  Tp  equal  to  3  sec  (see  Eq.  4)  is  used 
as  f  (t)  . 

In  Fig.  7-a  Is  shown  the  synthesized  seismic  image  produced  by 
the  vertically  Incident  plane  S  wave.  The  part  of  the  Image  produced  by 
the  direct  rays  gives  straightforward  Indications  about  the  structure. 
Suppose  this  image  is  a  recorded  one  instead  of  a  synthesized  one;  the 
following  procedure  is  suggested  to  invert  the  structure. 

First,  construct  a  primary  structure  model  based  on  the  part  of  the 
Image  produced  bv  direct  rays,  and  make  a  synthesized  Image  for  this 
model  to  compare  with  the  observed  image.  The  multiple  bounces  then 
provide  the  clues  for  refining  or  modifying  the  current  model.  This 
procedure  is  repeated  until  consistency  is  obtained. 

ii)  Example  2 

The  only  difference  between  this  example  and  the  last  example  is 
the  incident  direction  of  the  plane  S  wave,  which  is  now. 


The  seismic  image  is  plotted  in  Fig.  7-b.  Because  the  incident 
wave  is  slanting,  the  image  of  the  structure  is  tilted  as  expected. 

The  characteristics  of  this  image  actually  is  very  similar  to  that  of 
the  last  example.  For  instance,  the  triplication  feature  of  the  second 
multiple  in  each  case  appears  in  a  similar  way. 

iii)  Example  3 

The  only  difference  between  this  example  and  the  last  example  is 
the  geometry  of  the  interface.  In  this  example,  the  representation  of 
the  Interface  is  of  the  form  of  (F.q.  2)  with  the  following  parameters. 


d  km,  c  =  5  km.  w  «  50  km. 


Notice  in  this  case  that  the  layer  is  thin  in  the  central  part  and 
becomes  thick  and  flat  at  the  edges. 

The  seismic  image  is  plotted  in  Fig.  7-c.  By  comparing  Fig.  7-b 
and  7-c,  we  can  see  the  images,  especially  the  parts  produced  by  mul¬ 
tiples,  are  strongly  revealing  the  difference  between  the  two  struc¬ 
tures.  Their  characteristic  patterns  are  certainly  clues  for  identi¬ 
fying  the  features  of  the  subsurficial  geometry. 


VI.  DISCUSSION  AND  CONCLUSIONS 


In  a  receiver  structure  involving  curved  boundaries,  the  wave¬ 
form  complications  are  mainly  due  to  the  concentration  and  triplica¬ 
tion  of  the  multiple  bounces.  From  the  case  of  a  one-dimensional 
structure  to  the  case  of  a  two-dimensional  structure,  the  variation 
of  the  wave  forms  is  very  apparent.  However,  from  the  already  com¬ 
plicated  case  of  a  two-dimensional  structure  to  the  case  of  a  three- 
dimensional  structure,  the  further  complications  are  existing  delicately 
in  the  coda,  as  shown  in  Case  1  and  2  of  Section  IV.  As  we  can  see,  it 
is  rather  hard  to  make  a  starting  model  for  inverting  a  3-dimensional 
structure  if  there  is  only  a  single  station  in  the  receiver  region. 

To  practice  the  teleseismic  prospecting  efficiently,  we  need 
seismic  arrays.  As  shown  in  Section  V,  we  have  developed  a  scheme 
to  synthesize  the  seismic  images,  in  which  the  direct  rays  associated 
with  multiple  bounces  can  strongly  characterize  the  subsurficial 
features  of  the  receiver  structures.  The  scanning  procedure  costs 
Just  about  the  same  in  the  case  of  a  single  station  and  in  the  case 
of  one  hundred  stations.  Of  course,  the  fine  tuning  procedures  cer¬ 
tainly  make  the  latter  more  expensive. 

On  the  other  hand,  the  accuracy  is  also  depending  on  the  fine 
tuning  procedures.  However,  this  scheme  is  designed  to  allow  the 
sacrifice  of  accuracy  in  exchange  for  the  computing  cost.  Actually, 
the  wave  forms  shown  in  Section  V  were  synthesized  with  less  accuracy 
than  those  in  Section  IV. 

There  are  some  other  inaccuracies  related  to  the  caustics  and 
the  critical  and  post-critical  reflections,  since  the  glorified 
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optics  Is  basically  a  high-frequency  approximation.  By  comparing 
with  the  results  of  the  finite  element  method,  the  accuracy  Is 
still  satisfactory  even  when  the  dominant  period  is  as  long  as  20 
sec  (see  Kong  and  Helmberger,  1978)  . 

With  the  development  of  the  glorified  optics  and  the  scheme 
for  synthesizing  seismic  Images,  with  the  plentiful  earth  quakes 
and  nuclear  events,  and  with  the  ever  expanding  and  improving  seismic 
array  techniques,  we  believe  the  teleseismlc  prospecting  is  a 
promising  tackle  to  obtain  the  precise  features  of  structures  like 
sedimentary  basins,  salt  domes,  geosyncllnes ,  continental  shelves, 
subduction  plates,  anomalies  on  Moho,  etc. 


\ 
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FIGURE 

Fig.  1 


Fig.  2 

Fig.  3 


Fig.  4 


CAPTIONS 

Each  of  the  four  cases  in  this  figure  shows  the  projections 
of  the  paths  of  multiple  bounces  which  are  started  along  a 
baseline  with  a  given  azimuth  in  the  scanning  frame.  The 
left  column  shows  the  projections  on  the  plane  y  *  0.  The 
right  column  shows  the  projections  on  the  plane  z  -  0.  The 
dashed  curves  in  the  right  column  represent  the  contours  of 
the  interfaces. 

This  figure  shows  the  horizontal  projections  of  eight  groups 
of  rays  in  a  scanning  procedure.  Each  group  of  rays  is 
started  along  a  baseline  with  a  given  azimuth  in  the  scanning 
f  rame . 

This  figure  shows  the  projections  of  ray  tubes  in  a  two- 
dimensional  structure  in  (a)  and  (b),  respectively.  The 
left  column  shows  the  projections  on  the  plane  y  «  0.  The 
right  column  shows  the  projections  on  the  plane  z  -  0.  The 
dashed  curves  in  the  right  column  represent  the  contours  of 
the  interfaces. 

This  figure  shows  the  geometry  of  a  structure  of  the  type  of 
a  two-dimensional  basin  over  a  half-space.  The  left  plot  is 
its  cross  section  on  the  plane  y  =  0.  The  shorter  dashed 
curves  in  the  right  plot  represent  the  contours  of  the  inter¬ 


face,  which  is  infinitely  long  in  the  v-direction. 
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Fig. 


Fig.  6 


Fig.  7 


This  figure  shows  the  displacement  wave  form9  for  the  cases  in 
Section  IV.  The  left  column  shows  the  longer  period  response. 
The  right  column  shows  the  shorter  period  response.  The 
letters  x,  v,  z  are  indicating  which  components  they  are. 

This  figure  shows  the  geometry  of  a  structure  of  the  type 
of  a  three-dimensional  basin  over  a  half-space.  The  left  plot 
is  its  cross  section  on  the  plane  y  =  0.  The  shorter  dashed 
curves  in  the  right  plot  represent  the  contours  of  the  inter¬ 
face  . 

The  synthesized  seismic  image  for  the  examples  In  Section  V. 
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The  arriving  time  (T),  the  spreading  effect  (S) ,  and  the 
reflect  Ion- transml salon  effect  (J)  of  each  arrival  In 
Case  1  of  Section  IV.  The  number  of  segments  along  the 
geometric  path  of  an  arrival  is  Indicated  in  the  second 
column.  The  mode  of  each  arrival  is  also  indicated  there. 
Notice,  the  P-S  converted  phases  are  not  considered  here. 

The  arriving  time  (T) ,  the  spreading  effect  (,S) ,  and  the 
reflection-transmission  effect  CJ)  of  each  arrival  in 
Case  2  of  Section  IV. 

The  arriving  time  (T) ,  the  spreading  effect  (S) ,  and  the 
reflection-transmission  effect  CJ)  of  each  arrival  in 
Case  3  of  Section  IV. 

The  arriving  time  (T) ,  the  spreading  effect  (,S) ,  and  the 
reflection-transmission  effect  CJ)  of  each  arrival  in 
Case  4  of  Section  IV. 

The  arriving  time  (T) ,  the  spreading  effect  (S) ,  and  the 
reflection  transmission  effect  CJ)  of  each  arrival  in 
Case  5  of  Section  IV. 
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IV.  SLIPPING  INTERFACES:  A  POSSIBLE  SOURCE 
OF  S  RADIATION  FROM  EXPLOSIVE  SOURCES 

By  Carlos  Salvado  and  Bernard  Minster 

Introduction: 

Considerable  effort  has  been  devoted  to  the  understanding  of  SH 
radiation  from  underground  explosive  sources.  A  variety  of  possible  mechanisms 
have  been  proposed,  as  reviewed  by  Akl  and  Tsai  (1972).  They  Include  (1) 
mode  conversion  at  Irregular  Interfaces  (e.g.,  Akl  and  Lamer,  1970),  (2) 
cracking  In  the  vicinity  of  the  shot  point  (e.g.,  Klsslinger  ejt  al . ,  1961), 

(3)  prestress  relaxation  associated  with  the  creation  of  the  cavity  (e.g., 
Archambeau,  1972,  1973)  and  (4)  triggering  of  an  actual  earthquake  with  large 
scale  faulting,  a  model  favored  by  Akl  and  Tsai  (1972). 

It  Is  likely  that  a  combination  of  these  mechanisms  acts  to  produce  the 
observed  SH  and  Love  wave  radiation  from  underground  explosions.  The  relative 
contributions  of  various  mechanisms  will  of  course  depend  on  local  characteristics 
of  the  site.  Nevertheless,  It  may  be  noted  that  the  two  last  mechanisms  listed 
above  involve  the  release  of  strain  energy  stored  In  the  medium  prior  to  the 
experiment,  and  therefore  require  the  presence  of  sufficient  prestress.  On  the 
other  hand,  the  two  first  mechanisms  can  operate  even  in  the  absence  of 
significant  Initial  stress,  as  evidenced  by  the  observations  of  Klsslingei  et  al. 
(1961).  These  authors  observed  SH  radiation  generated  by  small  scale  explosions 
detonated  In  mud,  and  suggested  that  near  source  cracking  was  the  most  likely 
cause  of  shear  wave  radiation. 

In  order  to  better  understand  the  Influence  of  site  characteristics  and 
prestress  on  SH  and  Love  wave  generation,  an  Improved  assessment  of  the  phenomena 
which  do  not  require  prestress  is  needed.  Cracking  and  block  motions  in  the 
vicinity  of  explosions  fall  in  this  category;  a  review  of  these  problems 
Is  provided  by  Bache  and  Lambert  (1976). 

The  purpose  of  this  paper  is  to  demonstrate  that  the  presence  of  impeifectly 


ii 


bonded  interfaces  and  joints  near  a  purely  dilational  source  could  account  for 
far  field  transverse  components  of  displacement,  and  thus  with  a  suitable 
geometry,  for  SH  and  Love  wave  radiation  (e.g.,  Salvado  and  Minster,  1977). 

The  problem  is  greatly  simplified  because  we  restrict  ourselves  to  linear  boundary 
conditions,  and  solve  the  problem  in  the  far-field,  first-motion  approximation. 

It  appears  that  if  the  bonding  is  sufficiently  weak,  and  if  the  geometry  is 
favorable,  at  least  a  portion  of  the  observed  shear  radiation  could  be 
explained  by  this  mechanism. 


1. 


1.  Statement  of  the  Problem,  and  Boundary  Conditions 

Consider  the  simple  geometry  depicted  on  Figure  1:  ■  purely  dllatatlonal 

point  source  Is  located  at  a  height  h  above  the  plane  Interface  between  two 
elastic  half  spaces.  To  simplify  the  problem,  we  restrict  our  attention  to 
the  case  when  the  two  half  spaces  are  identical.  The  boundary  conditions 
usually  adopted  in  seismology  Involve  continuity  of  tractions  and  dlnolacements , 
and  In  this  case,  the  problem  reduces  to  that  of  a  point  source  In  an  Infinite 
space. 

We  shall  now  relax  the  boundary  conditions,  and  request  continuity  of  tractions 
and  of  normal  displacement  but  allow  a  jump  in  the  tangentlonal  displacement 
In  other  words,  slippage  is  allowed  between  the  two  half  spaces.  Sezawa  and  Kanal 
(1940),  Kanal  (1961)  and  Murty  (1975,  1976)  proposed  a  boundary  condition  which  does 
allow  such  slippage  to  occur,  and  has  the  advantage  of  being  linear;  it  is 
best  described  as  the  result  of  the  following  limiting  process. 

Suppose  the  two  half  spaces  are  actually  coupled  through  a  viscous  layer 
of  thickness  H  and  viscosity  n.  The  geometry  is  described  on  Figure  2.  Assuming 
a  plane  Couette  flow  In  the  viscous  layer,  the  shear  traction  is  related  to  the 
jump  in  tangential  velocity  by 

ori.  n^»j|ir(r,0+,t)  -ir(r,0',t)|  (0  <  n  <  ->  (1.1} 

We  now  allow  the  thickness  H  and  the  viscosity  n  to  go  to  zero,  but  require  that 
their  ratio  remain  constant.  For  convenience  in  the  subsequent  analyses,  we 
define 


a-  ♦ ...  a 
h  i  -  *  e 


(O' <  *  <  1)  (1-2) 


whore  u  is  the  rigidity  and  B  the  shear  wave  velocity  in  the  two  half  spaces. 
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♦  is  a  dimensionless  bonding  parameter.  In  this  limit,  the  boundary  condition 
becomes 


whe  re 


Nj 


i  I  V1"*0,0 


is  the  Jump  of  tangential  velocity  across  the  interface: 


|  llr  (r,0,  t)  +  =  iir(r,ott)  -  iir(r,0,t) 

When  ♦  ■  0,  the  tangential  raction  vanishes  and  the  interface  is  perfectly 
lubricated;  when  ♦  *  1,  no  relative  motion  between  the  two  half  spaces  is 
allowed  and  the  interface  is  perfectly  bonded. 

The  main  advantage  of  the  boundary  conditions  (1-3)  is  its  linearity, 
which  greatly  simplifies  wave  propagation  problems.  Imperfectly  bonded  inter¬ 
faces  have  been  modelled  using  nonlinear  boundary  conditions  (e.g., 

Chez  et  al.  (1978)  and  Miller  (1979)),  but  the  difficulties  are  such 
that  only  plane  waves  have  apparently  been  considered. 


3. 


and  simple  closed  form  solutions  may  be  found  in  the  first  motion  approximation. 
We  consider  the  potential  decomposition  (e.g.  ,  Cagniard.  1962) 
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and  the  solutions  of  the  form 
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where,  for  t  >  0,  the  source  medium, 


and.  for  z  <  0: 
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and  a,  0  are  the  P  and  S-vave  velocities  respectively. 

The  unknown  coefficients,  (Q1(k),  Sj(k);  i  ■  1,2)  are  to  be  determined  by 
use  of  the  boundary  conditions.  Continuity  of  tractions,  normal  displacement. 


slip  condition  (1-3)  lead  to  the  following  system  of  simultaneous  equations. 
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The  determinant  In  (2-5)  is  given  by 
A(k)  •  4Pvn  a2  D(k) 


D(k) 
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Notice  that  the  bracketed  expression  in  D(k)  is  the  Rayleigh  denominator 

The  solutions  to  (2-5)  are 

-v„h 
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In  the  case  of  a  perfectly  bonded  Interface,  ♦  ■  1,  and  the  only  nonvaniahing 
coefficient  is  Q.,(k),  as  expected. 

By  considering  derivatives  of  theae  coefficients  with  respect  to  the  bonding 
parameter  *,  one  finds  that  no  extrema  occur  for  *  ]0,1[.  In  other  words 
the  most  efficient  conversion  of  P  waves  to  S  waves  must  occur  for  ♦  ■  0,  that 
is,  a  perfectly  lubricated  interface. 

We  can  now  substitute  these  coefficients  in  the  solutions  (2-3)  and(2-A) ,  and 
replace  k  in  terms  of  the  ray  parameter  p  (e.g.,  de  Hoop,  1960;  Helmberger, 

197A;  Langston,  1976). 

k  -  -i  s  p.  (2-9) 


After  application  of  the  Schwartz  reflection  principle,  the  solution  takes 


the  form 


where 
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and  where  Che  following  definitions  hold 
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P(p)  -  l(p  -  2 up2)2  +  4u2n  URp2)  (1-*)  +  ♦  (2-13) 

a  u  p. 

Asymptotic  approximations  of  (2-10)  and  first  motion  approximations  can 
then  be  obtained  in  the  usual  fashion  (e.g.,  Langston,  1976).  By  substitution 
of  the  asymptotic  results 


K  (spr)  ~ 


ti  c~spr 
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the  solution  (2-10)  are  reduced  to  integrals  of  the  form 


l(r,z,s)  *  Je~st^r  •"t,1C^f(p)dp 


(2-15) 


from  which  the  displacement  must  be  retrieved  via  the  formulae  (2-1).  This 
requires  finding  the  spatial  derivatives  of  I.  Keeping  only  the  far-field 
terms,  we  have 
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3r  J  xsr  J 


r‘- 

*m,z**2  ,  ,a  /X  /e~8t(r,l,p)  |”f(p)dp 

3z  y-  y 


C2-17) 


The  body  wave  contributions,  which  we  are  Interested  in,  are  obtained  by  changing 
the  variable  of  integration  to  t  so  that  by  oae  of  Cauchys-Goursats  theorem 


on  the  complex  t-plane  we  have 
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3r 
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The  geometrical  arrival  time  ia  defined  by 


lA  *  t(r.*.PD> 


ll  -0 

3P  P0 


(2-20) 


If  F(t)  is  the  source  function  for  the  potentials,  then  the  time  domain 
equivalents  to  (2-18)  and  (2-19)  are 


Tt  r<t)  *  7^2  *  *"  (2-21) 

^  -  vV?  ft  f'(t)  '  fl72  “  '*  (2-22) 


where  *  indicates  a  convolution.  To  compute  the  first  motion  approximation 


we  note  that  t(r,z,p)  is  of  the  form 


t  -  pr  +  nagi(*)  +  n6g2(*)  (2-23) 

where  g1(i)  and  g2(r)  are  listed  in  Table  1  for  the  various  rays.  Using  the 
definition  (2-20)  and  replacing  t  by  tA  in  the  first  motion  approximation,  we 


;!£„!  2(t-t  )  «1(1> 


(2-24) 


so  that,  by  substitution  in  (2-21)  and  (2-22),  and  reduction 
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The  displacement  components  for  the  various  rays  can  be  obtained  by  using  these 
results  in  (2-1)  and  substituting  the  quantities  from  Table  1. 


Table 
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3.  Radiation  Patterns 

Suppose  the  displacements  found  in  the  previous  section  are  rotated 
to  a  coordinate  system  that  renders  the  slipping  interface  vertical. 

For  the  purpose  of  our  block  motion  study,  such  coordinates  are 
Earth  coordinates. 

For  a  take-off  angle  of  90°  in  Earth  coordinates  all  S  radiation 
reflected  from  and  transmitted  through  the  interface  is  pure  SH  relative 
to  the  surface  of  the  Earth.  Displacement  radiation  patterns  in  that 
case  are  simple  to  understand  and  were  computed  for  different  values  of 
the  bonding  parameter  and  the  following  values  of  the  remaining  parameters: 


Parameter  Value 

a/8  1.72 


P 

h 


3 

3.2  gms/cm 
1.0  km 


R«[r2+(z-h)2]1/2 


1000.0  km 


Amplitudes  are  measured  perpendicular  and  parallel  to  the  ray  Joining 
the  receiver  and  the  source.  This  gives  rise  respectively  to  the 
transverse  and  radial  components  of  displacement  as  seen  in  Fig.  3.  The 
radiation  patterns  for  a  unit  step  source  function  are  given  in  Fig.  4. 
where  the  solid  and  dotted  lines  are  the  transverse  and  radial 
components  respectively,  and  where,  positive  for  the  transverse  component 
is  defined  to  be  counter-clockwise.  The  horizontal  axis  represents 
the  edge  view  of  the  slipping  interface  and  the  units  on  this  axis 
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represent  the  ratio  of  the  amplitude  of  the  component  to  that  of  the 
incident  P  wave. 

These  patterns  are  seen  to  change  shape  with  bonding.  For  <f>  «  0.0 
the  angle  of  maximum  displacement  corresponds  to  angles  of  critical  reflections. 
Similar  lobed  patterns  were  found  by  Burridge  et  al.  (1964)  for 
reflections  from  a  free  surface. 

There  is  good  conversion  for  bonding  <p  =  0.0:  the  amplitude  of 
the  SH  wave  is  at  least  0.5  of  the  P  wave.  Conversion  gets  progressively 
smaller  as  the  bonding  is  varied  to  higher  values.  For  <(>  =  0.5  and  0.9, 
the  amplitude  of  the  SH  wave  is,  respectively,  .1  and  .01  of  the  P 
wave  amplitude.  This  is  due  to  the  fact  that  as  the  bonding  becomes 
stronger  the  jump  in  shear  displacement  across  the  interface  decreases, 
and  vanishes  when  the  bonding  is  perfect,  i.e.,  0  =  1.0.  This  is  seen 
in  Fig.  5. 

So  it  appears  that  a  low  value  of  the  bonding  parameter  is  needed 

for  this  mechanism  to  be  an  efficient  shear  wave  generator.  However 

these  radiation  patterns  are  subject  to  change  in  shape  and  amplitude 

with  changes  in  the  values  of  the  parameters  in  the  expressions  for 

the  displacement.  Most  important  of  these  is  the  perpendicular  distance 

h  from  the  source  to  the  slipping  interface.  As  can  be  appreciated  in  the 

-1/2 

expressions  for  the  displacements,  2-25  and  2-26,  these  vary  as  h  in 

the  high  frequency,  first  motion  approximation.  Thus  this  linear 
theory  predicts  that  if  a  slipping  joint  is  sufficiently  close 
to  an  underground  explosion  it  could  be  an  efficient  converter  of  P  to 
S  waves.  One  must  remember,  however,  that  the  approximation  developed 
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in  the  previous  section  is  far  field  and  linear  and  therefore  must 
break  down  for  small  h,  so  that  this  conclusion  is  only  qualitative. 

This  faulting  process  due  to  the  presence  of  a  joint  in  the  vicinity 
of  an  underground  explosion  gives  rise  to  a  different  displacement  pattern 
than  that  given  either  by  the  relief  of  prestress  along  a  fault  or  the  impinging 
of  P  waves  on  a  slipless  impedance  mismatching  surface,  both  of  the  same 
orientation  as  the  joint.  As  appreciated  in  Fig.  4  the  SH  displacement  radiation 
pattern  has  a  cylindrical  symmetry  with  respect  to  the  line  through  the 
source  and  perpendicular  to  the  interface:  it  is  symmetric  with  respect 
to  the  interface  but  antisymmetric  with  respect  to  the  line  perpendicular 
to  the  interface  through  the  source.  The  relief  of  prestress  along  a  vertical 
fault,  in  this  case,  a  vertical  strike  slip  fault,  has  a  double-couple 
radiation  pattern  which  is  symmetric  with  respect  to  the  interface  and  the 
line  perpendicular  to  the  interface  through  the  source  as  can  be  appreciated 
in  Heaton  (1979).  On  the  other  hand,  the  radiation  pattern  from  the  incidence 
of  P  waves  on  a  slipless  impedance  mismatching  surface  has  also  cylindrical 
symmetry  with  respect  to  the  line  perpendicular  to  the  interface  and 
through  the  source  but  it  differs  from  the  radiation  pattern  of  the 
slipping  joint  in  that  it  is  in  general  asymmetric  with  respect  to  the  inter¬ 
face  as  is  seen  in  Ewing  et  al.  (1957). 

This  distinction  of  radiation  patterns  may  provide  the  best  tool 
through  measurement  to  determine  whether  block  motion  is  indeed  a 
significant  contributor  of  S  waves  from  large  underground  explosions. 

It  does  not  matter  that  the  actual  rheology  of  the  fault  is  non-linear, 
for  example  Coulomb  friction,  because  the  radiation  pattern  of  the  first 
motion  will  have  the  same  polarity  as  in  this  study.  On  the  other  hand 
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If  the  test  site  Is  significantly  prestressed,  observation  of  block  notion 
through  the  polarity  of  the  radiation  pattern  may  be  difficult  to 
observe.  Since  Novaya  Zemlya  exhibits  no  seismicity,  and  this  points 
to  small.  If  any  prestress,  the  records  from  underground  explosions  at 
that  site  would  be  those  Indicated  to  test  this  theory. 

4.  _  Summary  and  Conclusion 

We  have  shown  that  a  slipping  Interface  separating  Identical  materials 
give  rise  to  conversion  of  P  to  S  waves.  This  conversion  is  best  when  the 
materials  are  perfectly  uncoupled  and  when  the  distance  from  the  source 
to  the  Interface  Is  small.  For  these  reasons,  block  motion  may  be  an 
efficient  generator  of  S  waves  from  large  underground  explosions  If  cracks 
are  quite  near  the  source  and/or  the  friction  on  the  crack  faces  is  small.  To 
be  able  to  determine  whether  this  model  Is  valid,  one  can  make  use  of  the 
displacement  radiation  patterns  from  tests  in  nuclear  sites  that 
exhibit  no  seismicity,  as  Novaya  Zemlya. 
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Figure  Captions 


Figure  1.  A  purely  dllatational  point  source  Is  located  at  a  height 
h  above  the  slipping  interface  separating  Identical  half¬ 
spaces. 


Figure  2.  The  Identical  half-spaces  are  separated  by  an  interface  which 
is  assumed  a  fluid  filled  gap  of  viscosity  n< 


Pigure  3.  Amplitudes  are  measured  perpendicular  (transverse  component 
T)  and  parallel  (radial  component  R)  to  the  ray  joining 
source  with  receiver. 


Figure  4.  Radiation  patterns  for  a  take-off  angle  of  90°.  Solid  line 
represents  the  transverse  component  while  the  dotted  line 
represents  the  radial  component. 


Figure  5.  Displacement  jump  across  the  slipping  interface  for 
various  values  of  bonding  parameter  $. 
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